Momentum-space interferometry with trapped ultracold atoms 
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Quantum interferometers are generally set so that phase differences between paths in coordinate 
space combine constructive or destructively. Indeed, the interfering paths can also meet in mo- 
mentum space leading to momentum-space fringes. We propose and analyze a method to produce 
interference in momentum space by phase-imprinting part of a trapped atomic cloud with a detuned 
laser. For one-particle wave functions analytical expressions are found for the fringe width and shift 
versus the phase imprinted. The effects of unsharpness or displacement of the phase jump are also 
studied, as well as many-body effects to determine the potential applicability of momentum-space 
interferometry. 

PACS numbers: 37.25.+k, 03.75.Dg, 42.50.-p 



I. INTRODUCTION 

In most quantum interferometers the phase differen- 
tial between paths that meet in a coordinate-space point 
or region at a given time lead to constructive or de- 
structive wave combinations and thus to fringes, but the 
paths can also interfere in momentum space and pro- 
duce momentum-space fringes. In particular, during the 
crossing of a wavepacket over a small and thin barrier 
(compared to the energy and width of the wavepacket), 
see Fig. QJ,, the momentum distribution can change dra- 
matically, vanishing at the center of the distribution, and 
being enhanced at the wings. This process would violate 
classical energy-conservation [l], and is due to interfer- 
ence in momentum space between incident and transmit- 
ted parts of the wave @, Q . 

The experimental implementation of this effect is in 
principle possible with current cold-atom technology, by 
turning off an effective detuned-laser barrier in the midst 
of the wavepacket passage, but a version which is simpler 
to implement is described here. The effect of the scatter- 
ing barrier in the original proposal is to imprint an ap- 
propriate phase on approximately half of the wavepacket, 
and this can also be achieved by shining part of a trapped, 
initially stationary, wavepacket with a strong laser pulse, 
during a short time in the scale in which a perturba- 
tion propagates (the correlation time Q), see Fig. QJ>. 
The phase imprinting technique was first introduced with 
the purpose of generating vortices Q . Here we shall 
study the properties of the resulting fringes, and show 
that the effect on the momentum distribution is similar 
to the effect of the scattering process, with the creation 
of a vanishing point ( "dark notch" ) at the center and en- 



(a) 



(b) 



9 




Trap 



FIG. 1: Schematic settings: (a) phase imprinting caused by 
wave packet passage above a weak, narrow potential, (b) 
phase imprinting caused by illumination with a detuned laser 
of a trapped wave packet. 
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hancement of the wings. Furthermore, we shall study the 
shift, width and visibility of this central "dark notch" as 
a function of the imprinted phase, as a necessary step 
to determine the potential applicability of momentum- 
space interferometry. The imprinted phase carries in- 
formation about the laser interaction (time, laser inten- 
sity, frequency) that can be obtained from the notch. 
By immediately removing the external trap, the momen- 
tum distribution is essentially frozen after the imprint- 
ing, and many-body effects cease to play a role. Then, 
the momentum-space notch will become by expansion a 
coordinate-space notch measurable with standard time- 
of-fiight techniques. Alternatively, the momentum distri- 
bution can be accessed by stimulated Raman transitions 
hj or through the single-particle reduced density matrix 



In the following section we will describe the setting in 
more detail for a single particle or many non-interacting 
particles. In Section Hill we will consider the role of in- 
teractions within the mean-field regime and in Section 
IIVI we will look at Tonks-Girardeau and non- interacting 
Fermi gases. 
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II. NON-INTERACTING REGIME 



which has a zero at qo, a solution of 



We start considering the non-interacting regime in 
which the single-particle description is valid. A highly 
anisotropic three-dimensional harmonic trap is assumed, 
so that the transverse degrees of freedom remain frozen 
and the system becomes effectively one-dimensional, 
along the axis with lowest trap frequency u>. It is useful 
to introduce dimcnsionless variables, namely a dimcn- 
sionless position y = y^pa;, where x is the dimensional 
position and m the mass of the single particle, and a 

dimensionless momentum q — w-^-fc, where k is the 
dimensional wavenumber. The Hamiltonian describing 
the system is H = —^—-§-^1 + ^^^x and in the above 

J 2m ax A 2 

dimensionless variables we get H = 4=r {— ^ 

The corresponding eigenvalues are E n = A=r(2ra 

the eigenstates are ip n (y) — , 1 H n (y)e~ y2 / 2 , 

H n (y) is the n-th Hermite polynomial. The eigenstates 
are normalized such that J dy ip n (y)ip m (y) = <5n,m- 

Initially, the trapped particle is described by the wave- 
function ipo(y) in coordinate space. Then a phase is im- 
printed on the right hand side of the trap, i.e. for y > 0. 
For atoms, this can be achieved by shining an appropriate 
detuned laser pulse for a short time t. The detuned laser 
acts as a mechanical potential V® (y) on the atom, where 
V = f2 2 ?i/4A, ft is the Rabi frequency, and A the detun- 
ing (laser frequency minus transition frequency) . If the 
time t is short, the effect is to imprint a phase tp = — Vt/% 
on the wave function for y > 0. The wavefunction in co- 
ordinate space becomes ^p (y)e lvw ^ with w(y) — 0(y), 
and in momentum space 
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dyMv)^ iqy 
dy ip (y)e~ iqv . 



Each momentum gets an amplitude contribution from 
two different terms and we may expect interferences in 
I0o(<z)| 2 for tp > 0. In the following, this interference 
pattern will be studied and analytical expressions will be 
found for the fringe shift, width and visibility versus the 
phase tp imprinted. The effects of unsharpness or spatial 
displacement of the phase jump are also studied. 
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Momentum distributions for this reference case after 
different phase imprintings are displayed in Fig. [2^. Note 
the optimality of ip = tt to produce a deep minimum, in 
fact a zero, exactly at the peak of the original distri- 
bution, qo = 0, and the enhancement at the wings. In 
the following, we will concentrate on this central "dark 
notch" . The "motion" of qo with ip can be seen in Fig. 
[2^,, and in more detail in Fig. [3^ (solid line), where 
qo is plotted versus ip. It shows a linear behavior in 
7r/2 < tp < 3tt/2 and slight deviations beyond that range. 
The width Aq = q + — q_ of the central interference dip 
is defined as the difference between the momentum q + 
of the maximum on the right-hand side of the minimum 
and the momentum q_ of the maximum on the left-hand 
side. Fig. [3Jd (solid line) shows the width versus the 
phase. Note that the width is always greater that \/27r 
(value of the thick dotted line) . 

Another important quantity is the visibility of the min- 
imum, which we define as 



min± |<?!>o(?±)l ^ 100(90)1' 
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This visibility is plotted in Fig. [3J: (solid line). From 
the calculations we can infer that the visibility limits the 
working range of the interferometer to 7r/2 < <j> < 37r/2, 
and is optimal around cj) = n. 



B. Analytical approximations for the reference case 

The goal now is to derive approximate analytical for- 
mulae describing the properties of the central "dark 
notch" as a function of the imprinted phase. From 

= o, we get the extreme points of |</>o(<?)| 2 as 
solutions of 
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A. Reference Case 

Let us first study the effect of imprinting a phase tp on 
the ground state of the harmonic trap, n — 0. In this 
case the momentum probability density becomes 



\Mq)\ 2 = 



cos 



(I) + sin (|) erfi 



0. 
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Note that if q is a solution of for tp = ir + Atp then 
— q is a solution for tp = n — Atp. 

One of the solutions of equation © fulfilling z\ = 
and describing the momentum of the minimum is ap- 
proximately given by 
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FIG. 2: Wave function in momentum space for different phase 
imprintings tp; (a) reference case, yo = 0, ( = 0, (b) effect of 
shifting, yo = 0.3, £ = 0, (c) effect of smoothing, yo = 0, 
£ = O.f , (d) effect of interaction, solution of the GPE: yo = 0, 
C = 0, 9 = 20. 



(based on a linearization of z\ = around y> w tt and 
<7 w 0). This describes a linear displacement of the mini- 
mum with tp. Fig. shows the exact momentum of the 
minimum qo (solid line) and qo (thick dotted line) versus 

tp. 

Now we shall obtain expressions for the (momentum of 
the) left maximum, and the right maximum q + . An 
approximate solution of z 2 — is 



f 

/2tt 
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FIG. 3: (a) Momentum of the minimum versus tp, (b) width 
of the minimum versus tp, (c) visibility of the minimum versus 
tp; in all cases n = 0. 

Analytical approximations: thick dotted lines: v, Aq, qo; dot- 
ted lines: v, Aq, qo (Thomas-Fermi with g = 20). 
Exact results: solid line: reference case, v, Aq, qo with 
yo = 0, £ = 0; boxes: effect of shifting, v, Aq, qo with 
yo = 0.3, ( — 0; triangles: effect of smoothing, v, Aq, qo 
with yo — 0,£ = O.f; circles: effect of interaction, v, Aq, qo 
based on the solution of the GPE with yo = 0, ( = 0, g — 20. 



and q « 0. Another approximate solution of z 2 = is 

27T — tp 



/2tt 



obtained by linearizing z 2 = around ^ « 27r and g«0. 
Thus we get for the the width of the interference dip 



A<? = q+ - q- 



q- = V2tt =: Aq. 



which follows from a linearization of z 2 = around tp w 



Figure [SJd compares the numerically calculated exact 
width Aq (solid line) with Aq (thick dotted line). 

We can also find a simple expression for the visibil- 
ity. From Eq. ([1]) and using also the approximations of 
the momentum of the two maxima and the minima and 
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retaining only the first order in ip we get 
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The final result is 
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also shown in Fig. [3J: (thick dotted line). It gives a lower 
bound for the exact result v (solid line). 



C. Perturbations of the reference case 

In this subsection we examine the effect of perturba- 
tions of the reference case. First we want to discuss the 
effect of shifting the edge of the phase imprinted region, 
yo, out of the center of the trap, i.e. we have w(z) — 
©(y — Ho)- Right after the phase-imprinting, the 0-th 
eigenstate now becomes ip (y)e i<pw W = ip (y)e kfle ^- y °\ 
and the momentum distribution becomes 
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dq 

Of 
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: 0, we get the extreme points as solutions 
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where z — cos (^) + sin (^) erfi ( 9 Jf° ) ■ Again, if q is 

a solution of © for ip = 7r + Aip then — q is a solution 
for ip = 7r — Aip. In addition, because ([6]) is not changing 
if yo is replaced by —yo, the solutions are the same for 
±yo. If yo <C 1, to first order in yo Eq. (O and Eq. ([6|) 
are independent of yo and therefore the above derived 
approximations in Section III Bl for the case yo = still 
hold. An example for yo — 0.3 is plotted in Fig. [2Jd. 
The corresponding momentum, width, and visibility of 
the minima versus ip is also plotted in Fig. [3] (boxes). 
The main effect of increasing yo is to lower the visibility 
(Fig. 0b). 

Finally, we want to look at the effect of a more realistic 
smooth profile of the imprinted phase, instead of using 



an idealized step function. Therefore, we consider now a 
sigmoid function 



w{y) = ~[l + tanh(y/0], 



(7) 



which for £ — > becomes O(y). The results for a smooth- 
ing £ = 0.1 can be seen in Fig[2j: and in Fig. [3] (triangles) . 
Smoothing results mainly in a shift of the maximum of 
the visibility (see Fig. [3];). 



D. Momentum interference for Excited States 

We shall next consider the effect of phase imprinting 
on excited states of the harmonic trap with the simplest 
profile w(y) = 6(y). The probability amplitude in mo- 
mentum space is then given by 



i 







27r2 n n!V7F 
dyH n (y)e-y 2 / 2 [(-1) 



which clearly simplifies for n = to Eq. (TT]). 

Let us look for q fulfilling \(f> n (<la)\ — 0, i.e. for the 
momentum of the minimum. Assuming q <C 1 such that 

(_l)n e iMo +e i¥> e --iOTo ~ ({-l) n +e l v)+iqo{{-l) n -e iv )y, 

1 
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*v 



)Bn\ 



where we have introduced A n = J Q °° dy H n {y)e y2 / 2 : and 

B n = f Q °° dyyH n (y)e~ y2/2 . Solving this for |0„(go)| 2 = 
0, we get 



. (-!)" + e*f A n 
\-l) n - e*v B n 



(8) 



The cases in which n is even or odd will be examined 
separately. 

a. n even: We are interested in the motion of the 
zero |0 n (go)| 2 = for ip « it. From Eq. ||5J) we get in 
first order in tp — n that 

A n ip-TT 



Qo 



B„ 



qo- 



Examples for the exact solution qo and the approximation 
qo for n = 0, 2 can be found in Figure 0Jl 

b. n odd: Now we are interested in the motion of 
the zero \(j> n {qo)\ 2 = versus ip for ip w 0. From Eq. §8§ 
we get in first order in ip that 

A n ip 

9o ~ 77-77 = : go- 



Examples for the exact solution q and the approxima- 
tion qo for n = 1, 3 can be found in Figure \Qp. 
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FIG. 4: Effect for excited states: (a,b) Momentum of min- 
imum versus tp; (a) n — 0: qo (solid line), q~o (thick dotted 
line); n = 2: qo (triangles), qo (dotted line); (b) n = 1: qo 
(solid line), qo (thick dotted line); n = 3: go (triangles), qo 
(dotted line); (c) Ratio versus n. 



In addition, Fig. 0b shows the value of the ratio 4 21 
for odd and even n. Clearly, increasing n makes the 
interferometer less sensitive to phase variations and n = 
provides the optimal behavior. 



III. THE MEAN-FIELD REGIME 

We shall consider now the role of the interactions 
within the mean-field approach. In the mean-field 
regime, for low enough temperatures the phase fluctua- 
tions can be suppressed jH . Weakly interacting ultracold 
gases in ID are then described by the Gross-Pitaevskii 
equation (GPE). 

Assume that an effectively ID Bose-Einstcin conden- 
sate is prepared in the harmonic trap. The condensate 
wave function is the ground state of the ID (stationary) 
GPE 



h 2 d 2 V(x) 
2m dx 2 
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FIG. 5: Interference in the momentum distribution of a Bose- 
Einstein condensate in the Thomas-Fermi regime. 



where /i is the chemical potential and g\u the effective 
ID coupling parameter related to the three-dimensional 
scattering length We assume that / dx |\l/(a;)| 2 = 

1. By introducing u :— 2fi/(?ud), g := \J~^giD, and 



ip(y) 



vnu) 



-^-y 1 we can write this equation in 



dimcnsionlcss form, 



ii>(y) = - 



d 2 ^(y) 
dy 2 



y 2 + g\i>(y)\ i>{y), 



with J dy \ip(y) | 2 = 1. The ground state can be nu- 
merically computed by using the imaginary time method 
[Til Il2l ]. For <p = 0, it is well known that as the mean- 
field interaction is increased, the density profile becomes 
more uniform, while the resulting momentum distribu- 
tion 1 0| 2 is sharply peaked [H,[3|- To study the effect 
of a small g as a perturbation of the previous results, we 
shall imprint a phase tp on the ground state wavefunction 
and calculate the minimum of the resulting interference 
pattern in momentum-space. An example with (7 = 20 is 
shown in Fig. The visibility, the width and the mo- 
mentum of the minimum versus ip with g = 20 is shown in 
Fig. [3] (circles). The main effect is that the slope in Fig. 

decreases with increasing atom-atom interaction, i.e. 
the sensitivity of the interferometer with respect to ip de- 
creases with increasing g. A large atom-atom interaction 
g may also perturb the measurement of the momentum 
distribution by time-of-flight techniques. There is how- 
ever also a positive effect: an increase of g makes the 
interference dip sharper and improves the visibility, see 
Fig. [3b. 

It is possible to derive analytical approximate formu- 
lae for large g. For g ^> 1 the condensate enters into the 
Thomas- Fermi regime [l3|, [l4| . The mean- field interac- 
tion is then so large that the kinetic energy can be ne- 
glected in the Hamiltonian so that the time-independent 

GPE reads utp(y) = (y 2 + g |^(y)| 2 ) ?/%)■ The 
Thomas- Fermi wavefunction is then given by ipTF{y) = 
y/(u — y 2 )/ g with u = (3g/4) 2 / 3 whenever \y\ < d and 
zero elsewhere; d — ^fu is the Thomas-Fermi half-width. 

The probability distribution in momentum space after 
a phase imprinting tp with a profile w(y) = 6(y) is given 
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in this case by 



\<hF(q)\' 



3vr 



J x {qd) cos (|) +Hi( 9 d)sin(|)] 2 .(9) 



Here, J\(k) is the Bessel function of first order and Hi(y) 
is the first order Struve function [l5|. \4>tf(q)\ 2 /d is plot- 
ted for different values of ip in Fig. [5] Again there is a 
minimum for ip — n at q = which is shifted if p is 
changed. 

The minima and the maxima of \4>tf{q)\ 2 for a fixed 
ip can be found by looking at the zeros of the derivative, 
this leads to the equation 



1 

qd 



Ji(<?d)cos ( ~) + Hi(gc2)sin (~ 



{qdf 



nqdJ2{qd) cos y— J + ^2 + TrgdH^^ci) sii 



= 0. 



Making a\ — and using a linearization around q ~ 
and <y9 « 7r, we arrive at 



3tt 
8d 



where d = (35/4) 1 / 3 , which allows to find d measuring 
the notch displacement. Making ai ~ and using a 
linearization around q « and ip « 0, we arrive at 



the similarities between the TG and spin-polarized non- 
interacting Fermi gases. The ground-state wavefunction 
of the later in a harmonic trap is the familiar Slater deter- 
minant, ® F (yi,...,y N ) = ■^jdet ( ^ k ^^ ) ip n {y k ) 1 built 
from the set of single-particle orthonormal eigenstates 
{ip n (y)}- Such atom Fock state can be efficiently pre- 
pared using the atom culling technique as described in 
[la ]. Note that the wavefunction typ is totally antisym- 
metric and vanishes whenever the positions of two parti- 
cles coincide. The TG wavefunction is obtained from W p 
by imposing the correct symmetry under permutation of 
particles, i.e., using the Fermi-Bose (FB) mapping [l6[ 

^Tc(yi, ■ ■ ■ ,Dn) = n sgn(y k -y j ) i &F(yi,...,y N ). 

l<j<k<N 

Clearly, both dual systems share the 
same density profile [17[ Ptg If (2/1 1) = 

N J \^TG/F{y, V2, ■ ■ • , yn] t)\ 2 dy 2 ■ ■ ■ dy N = 

J2n=o l^niUi t)\ 2 , a s it is the case for any other lo- 
cal correlation function. However, their momentum 
distributions 



n(q) = (2tt)- 



dydy e 



'pig(v-v') 



p(y,y') 



(10) 



are drastically different. Provided that the reduced 
single-particle density matrix (RSPDM) of spin-polarized 
fermions is 



PF{y,y') 



jV-l 

E 



fPnivWniy'), 



(ii) 



^^^^ 

and by using a linearization around q and p w 27r, 
we arrive at 



An estimate of the width is then Aq = q+ — q~- = || . An 
approximation for the visibility can be also derived as in 
Section MM 



v sa 1 — 0.5 \p — 7T I 



=: v. 



The approximate values of the notch momentum, width, 
and visibility in the Thomas-Fermi regime are also plot- 
ted in Fig. [3] (dotted lines). 



IV. THE TONKS-GIRARDEAU AND 
NON-INTERACTING FERMI GASES 

At low enough densities, and under tight-transverse 
confinement, ultracold gases enter the Tonks-Girardeau 
(TG) regime [16j . in which the strength of the effec- 
tive short-range interactions becomes so large that the 
mean- field theory fails [TtJ ■ Fortunately, Bose- Fermi du- 
ality offers a powerful and exact approach, exploiting 



the momentum distribution is the sum np(q) — 
J2n=o \^-n(q)\ 2 (where ij) n denotes the Fourier transform 
of ip n ). For the TG gas, an efficient way of computing 
the RSPDM has been introduced [H, 20], namely, 



N-l 

ptgM)= E Tp*(v) A in(y, v')i>n(y'), (12) 

l,n=0 

where A(y,y') = (P _1 ) T detP and the elements of 
the matrix P are Pi n (y,y') = J dzip*(z)4> n (z)sgn(z — 
y)sgn(z - y'), which reduces to Pi n = Si n - 

2 fjj dzipi(z)(j) n (z) for y < y' without loss of generality . 
The momentum distribution of the TG gas, can thus be 
obtained as a double Fourier transform. 

We consider a phase imprinting with w(y) = 6(y). Un- 
der this phase imprinting, a remarkable difference arises 
between the momentum distribution of both dual sys- 
tems. For a moderate N the visibility of the interference 
fringes in the TG gas is reduced (see Fig. [5^) but in 
the fermionic case the pattern has been washed out com- 
pletely (see Fig. [SJd). For larger N the visibility of the 
TG dip decreases. It is hence clear that the observa- 
tion of such effect in any of these dual systems would be 
difficult, and we turn our attention to a closely related 
alternative approach. 
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FIG. 6: Interference in momentum space, N = 10, a tp = tt 
phase is imprinted for y > (solid line), <p = case (dashed 
line); (a) TG gas, (b) Fermi gas, (c) TG gas after parity- 
selective evaporation, (d) Fermi gas after parity-selective 
evaporation. 



0.4 


i 

N=5 


1 1 


& 0.0 


_ N=10 

N=20 


» * * — - _ a 


-0.4 




1 



0.0 



0.5 



1.0 

<p[7t] 



1.5 



2.0 



FIG. 7: Displacement of the maximum of the momentum 
distribution as a function of the phase imprinted, for parity 
selective evaporation, fermionic cloud with only odd states. 
The symbols correspond to a smooth phase-imprinting profile 
with ( — 1/2, whereas for the lines £ = 0. 



Recently, a parity-selective evaporation (PSE) method 
has been proposed which allows to prepare in princi- 
ple excited states composed exclusively of odd-parity 
single-particle eigenstates [21] . This is achieved by 
shinning a blue-detuned laser at y ~ which re- 
moves the even-parity eigenstates. For a spin-polarized 
Fermi gas the excited many-body wavefunction becomes 



^f{Vi, ■ ■ ■ ,Vn) 



det^ fc=1 V>2n-i(2/fc)- The corre- 



/JV! 



sponding momentum distribution exhibits a well-defined 
zero at q = for all N which is stationary, and robust 
against significant smoothing of the phase-imprinting 
profile. The TG wavefunction equally follows from the 
Bose-Fermi map for PSE-prepared states. However, the 
n TG (<?) is qualitatively insensitive to the selected parity 
of the single-particle states, and lacks any principal peak 
(or dip) potentially useful for momentum-space interfer- 
ometry (see Fig. |BJ;). On the other hand, the pattern 
np(q) in the fermionic case is reversed under phase im- 
printing, turning a zero into a peak in the momentum 
distribution (see Fig. Let us consider again a phase 
imprinting with the sigmoid profile (J7J| . In Fig. [7] we 
have calculated the shift of the maximum in np(q) as a 



function of ip for the cases ( = (Heaviside function) and 
C = 0.5. The dependence is found to be linear even in the 
presence of the large smoothing in the profile C = 0.5. 

Therefore, between both dual systems the TG gas is 
preferred using phase imprinting, whereas in combination 
with PSE, the fermionic system is a better candidate. 



V. DISCUSSION 

The localized phase imprinting method 0] on trapped 
cold atoms has been discussed up to now mostly in con- 
nection with the generation and study of solitons. In 
this paper we have instead focused on the characteriza- 
tion of the momentum distribution right after the phase 
imprinting. 

First, phase imprinting of half of the wavepacket can 
be regarded as a simple way to realize the interferom- 
etry in momentum space that has been previously put 
forward for more complex scattering processes between 
cold atoms an weak laser barriers E| . Similar to the 
scattering setting, a central "dark notch" appears in the 
momentum distribution after phase imprinting, as well 
as an enhancement of the wings. An advantage with re- 
spect to the scattering method is that there is no need 
to make the width in momentum of the incident wave 
packet small to get the same transmission coefficient and 
therefore the same phase shift for all momenta. Thus 
we can make the trap tighter and tighter increasing the 
sensitivity. 

Furthermore, the characterization of the momentum 
distribution is a preliminary step to determine the po- 
tential applicability of momentum interferometry where 
an unknown phase should be determined from the mo- 
mentum shift of the central "dark notch" . 

We have studied different configurations, regimes, and 
perturbations. The momentum dark notch for non- 
interacting particles in the ground state provides the 
most sensitive meter for the imprinted phase among the 
different states considered. In dimensionless units, the 
momentum shift of the "dark notch" versus the imprinted 
phase is in this case approximately given by qo = \p^^-^~ 
(seeEq. (HJ). In dimensional units we get for the velocity 



vo 



flUl /7T^ 3 / 2 (fi — 7T 

~2, 



such that we can enhance the sensitivity concerning 
phase differences in principle to arbitrary high values by 
letting u> — > oo. If the external trap is immediately re- 
moved after the phase imprinting, the momentum dis- 
tribution is essentially frozen and the velocity vo can be 
measured with standard time-of-flight technique. Assum- 
ing a free time of flight of duration t after the phase im- 
printing the "dark notch" will move a distance s = tv . 
If we have a spacial resolution of As then the resolvance 
r of our momentum interferometer can be defined for the 



8 



reference case as 

7T t [hu> /7T\ 3 / 2 

r := ~Ktp = ~As\J ~m \2/ ' 

where Atp is the minimum resolvable deviation of the 
phase from n. For t = 200 ms, m =mass( 87 Rb), ui = 
2ir x 2 kHz and As = 5/nm, we get r ss 239. The effects 
of unsharpness or spatial displacement of the phase jump 
have also been studied and the results qualitatively still 
holds. Many-body effects in the mean-field regime lead to 
a mild sensitivity loss but also to an interesting increase 
of visibility. In all cases there is still a linear dependence 
of the "dark notch" velocity on the phase <p, i.e. vq 
a (f — 7r) such that the phase can be determined from 
the velocity of the "dark notch" . 

Other extreme regimes, as the Tonks-Girardeau gas 
of Bosons or an ideal Fermi gas diminish the interfer- 
ence, except, in the later case, when an auxiliary parity- 
selection procedure is applied to retain odd-states. A 
peak is then formed with linear dependence on the im- 
printed phase, very stable with respect to the smoothness 
of the profile of the imprinting laser. 



Finally, note that even though there arc no atom-atom 
interactions in the reference case, some of the aspects 
usually attributed to the solitons may already be recog- 
nized, in particular the formation of the dark notch in 
momentum representation and its shift with the value of 
the phase jump or its smoothness. Dynamical studies of 
the state evolution will provide further comparison with 
soliton dynamics and will be dealt with elsewhere. 
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